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I consider a stationary axisymmetric force-free degenerate magnetosphere of a rotating Kerr black 
hole surrounded by a thin Keplerian infinitely-conducting accretion disk. I focus on the closed-held 
geometry with a direct magnetic coupling between the disk and the event horizon. I first present 
a simple physical argument that shows how the black hole’s rotation limits the radial extent of the 
force-free link. I then conhrm this result by solving numerically the general-relativistic force-free 
Grad-Shafranov equation in the magnetosphere, using the regularity condition at the inner light 
cylinder to determine the poloidal current. I indeed hnd that force-free solutions exist only when 
the magnetic link between the hole and the disk has a limited extent on the disk surface. I chart 
out the maximum allowable size of this magnetically-connected part of the disk as a function of the 
black hole spin. I also compute the angular momentum and energy transfer between the hole and 
the disk that takes place via the direct magnetic link. I hnd that both of these quantities grow 
rapidly and that their deposition becomes highly concentrated near the inner edge of the disk as 
the black hole spin is increased. 


I. INTRODUCTION 


This paper is a short version of the full article M 

Magnetic fields play a crucial role in both the dy¬ 
namics and the energetics of accretion flows onto black 
holes, which makes it important to understand the 
global structure of an accreting black hole magneto¬ 
sphere. 

Conceptually, one can think of two basic types of 
magnetic-field geometry. The first type is the open- 
field conhguration responsible for the Blandford- 
Znajek mechanism for powering AGN jets 0. In this 
configuration there is no direct magnetic link between 
the hole and the disk: all the held lines are open and 
extend from the black hole and from the disk to inhn- 
ity. The open magnetic held extracts the rotational 
energy and angular momentum from a spinning black 
hole and transports them outward to power a jet. A 
similar process also works for the held lines connected 
to the disk. 

Recently, there’s been a rising interest in a different, 
namely, closed-held conhguration (Figure QJ, charac¬ 
terized by a direct m^netic coupliM between the hole 
and the disk iilliilllillllllllllll 
Here, energy and angular momentum are directly ex¬ 
changed between the hole and the disk via the mag¬ 
netic held. This exchange governs the spin evolution 
of the black hole H, [l3 and also has a strong effect 
on the energy dissipation prohle in the disk, which 
leads to important observational consequences 0 . 
In order to understand any of these processes, how¬ 
ever, one needs to have a good picture of the struc¬ 
ture of the magnetic held. Whereas there have been 
some numerical studies of open-held force-free con- 
hgurations Ellllllllll fl^ . the structure of the 
closed-held magnetosphere have remained essentially 
unexplored, with the exception of my recent work on 
the Schwarzschild case This provides the moti¬ 
vation for my present study p^. 



FIG. 1: Schematic drawing of a fully-closed black hole - 
disk magnetosphere. Energy and angular momentum are 
exchanged between the hole and the disk through a direct 
magnetic link. 


II. THE MODEL 


The main goal of this work is to determine the 
structure of the coupled disk-black hole magneto¬ 
sphere under the following assumptions: 1) everything 
is axially-symmetric and stationary; 2) the magneto¬ 
sphere is force-free and infinitely-conducting; 3) the 
disk is infinitely conducting, thin, and Keplerian. 

The structure of such a magnetosphere is governed 
by the generally-relativistic force-free Grad-Shafranov 
(e.g., 0, [ 23 ) for the poloidal flux function T (in 
Boyer-Lindquist coordinates): 
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Here, a is the lapse function, w is the cylindrical ra¬ 
dius, Di:’('I') is the held-line angular velocity, AD = 
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Dp — oj, where uj is the frame-dragging frequency, and 
I = /('f') is the poloidal current. 

This equation is very non-trivial, not only because 
it’s a non-linear partial differential equation, but also 
because it involves two, in general nonlinear, functions 
of 'h, namely and /('I'), that somehow need to 

be specified. 

In our particular problem, the first of these func¬ 
tions, is determined fairly straight-forwardly 

from the boundary conditions. Indeed, all the held 
lines threading the hole are also attached to an in- 
hnitely conducting Keplerian disk. Therefore, in 
steady state, the angular velocity of any given held 
line 'k is equal to the Keplerian angular velocity of 
this line’s footpoint ro('I') on the surface of the disk: 
Dp{^) = ('!')]. Thus, if a Dirichlet boundary 

condition for equation m on the disk surface is spec- 
ihed as 'k = d>d(r), one can invert '^d{r) and obtain 
an explicit expression for ^^’('I'). 

The poloidal current /(d>), on the other hand, can¬ 
not be determined as easily, since there is no surface 
on which one can just explicitly prescribe /('k). In¬ 
stead, one has to use an indirect procedure that makes 
use of the regularity condition at the inner light cylin¬ 
der. Indeed, as has hrst been noticed by Znajek M 
and by Blandford & Znajek Q, any rotating held line 
going into the hole has to pass through an inner Light 
Cylinder (LC), just as any rotating held line going 
out to inhnity has to pass through an outer LC. The 
light cylinder is dehned as the surface where the ro¬ 
tational velocity of held lines, as measured by a local 
Zero-Angular-Momentum Observer (ZAMO), is equal 
to the speed of light and where E = Bpo\ in the ZAMO 
frame: 


convenient for treating the LC regularity condition; in 
general, the LC sweeps through the grid until the hnal 
solution is achieved and so one has to actually evolve 
rLc(^) dynamically during the iteration procedure. 

Now let me discuss the role of the black hole event 
horizon, a = 0. The event horizon is also a regular 
singular surface of the Grad-Shafranov equation O- 
This is a very important observation; in particular, it 
means that one cannot, and indeed need not, provide 
any additional independent boundary condition here. 
Instead, one should impose the condition of regularity 
of the solution at the horizon. This condition can be 
actually obtained by simply setting a = 0 in equa¬ 
tion o, with the following result, commonly known 
as Znaiek’s event-horizon boundary condition (see 
also Q): 


/[vko(0)] 


2MrHSmO . ^ ddi'o 
_ tl _AO_ 1. 

’ 


r = rH, (3) 


where 'ko(^) = 'I'(r = rH,0) is the horizon magnetic 
flux distribution. In fact, this condition is what the 
fast magnetosonic critical condition becomes in the 
limit when the plasma density is taken to zero and the 
fast magnetosonic surface approaches the horizon m 
1^1^ . The signihcance of condition © is that it can 
be used to determine the magnetic flux distribution 
d>o(0) at the horizon provided that the functions /(d>) 
and n('k) are already known. 


III. MAIN EFFECT: DISRUPTION OF THE 
HOLE-DISK COUPLING BY THE BLACK 
HOLE’S ROTATION 


a = aLC = |AII|tx7 , (2) 

(in our problem, this inner LC is not at all a cylinder 
but rather an almost spherical surface between the 
event horizon and the boundary of the ergosphere). 
As one can easily see, the LC is a regular singular 
surface of equation ©■ This means that, in general, 
this equation admits solutions that are not continuous 
or continuously differentiable at r = rpc{0). That is, 
the regularity of the solution at the LC does not fol¬ 
low from the equation automatically, but has to be 
imposed as an additional physical condition. And it 
is precisely this condition that one can use to uniquely 
fix the function /(d>)! Actually doing this numerically 
requires an iterative approach in which /('k) is evolved 
together with (r, 9) until a regular solution of equa¬ 
tion o is obtained. This is very similar in spirit to 
the procedure used by Contopoulos et al. in their 
study of a force-free magnetosphere of an axisymmet- 
ric pulsar (see also mil ). One essential compli¬ 
cation in the black hole case is that the held lines 
rotate differentially and hence the actual location and 
shape of the LC are not known a priori. Therefore, 
one cannot set up the computational grid in a manner 


Another reason the event horizon regularity condi¬ 
tion is very important is that it enables us to arrive 
at a very important effect limiting the maximal possi¬ 
ble spatial extent of a direct force-free magnetic link 
between a rotating black hole and a conducting Kep¬ 
lerian disk. Here I present the basic physical reason 
for this effect (conhrmed by the numerical simulations 
presented in the next section). 

First, let us suppose that a force-free conhguration 
of Figure m where all the held lines attached to the 
disk at arbitrarily large radii thread the event hori¬ 
zon, does indeed exist. Then let us consider the polar 
region of the black hole, r = rp, 6 ^ 0, and suppose 
that near the rotation axis the hux 'ko(^) behaves as 
a power law: 'ko ~ (the most natural choice be¬ 
ing 'ko ~ Next, note that the held lines thread¬ 
ing the horizon in this polar region connect to the 
disk at some very large radii, ro(d>) ^ rp- For suf- 
hciently small and hence sufficiently large ro(d>), 
the Keplerian angular velocity is small, and therefore 
f2i;’('k) is much smaller than the black hole rotation 
rate Dp = aj^rp, so that AH —Dp = const. Then, 
from the event horizon regularity condition ||2Jl it fol¬ 
lows that /(d>) ^ —Dh'^ —a'k, as ^ 0, and 
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so 

^ 0) ~ (4) 

Now, if we look at the force-free balance on the same 
field lines but far away from the black hole, at dis¬ 
tances comparable with the footpoint radius, we can 
estimate the linear term on the left-hand-side of the 
Grad-Shafranov simply as We see that both 

the left-hand side and the right-hand side given by 
equation o scale linearly with dt but the left-hand 
side has an additional factor ^ r~^. Thus, at suf¬ 
ficiently large distances this term becomes negligible 
when compared with the J/'(4')-term given by 0). In 
other words, the toroidal field produced in the polar 
region of the horizon by the black hole dragging the 
field lines is too large to be confined by the poloidal 
field tension at large enough distances. In fact, this ar¬ 
gument suggests that the maximal radial extent i?max 
of the region on the disk connected to the polar region 
of a Kerr black hole should scale as i?max ~ rn/a in 
the limit a 0. In the Schwarzschild limit a —*■ 0, 
this maximal distance goes to infinity and hence a 
fully-closed force-free configuration can exist at arbi¬ 
trarily large distances, in agreement with my previous 
conclusions 0 . 

To sum up, even though the field lines can, to a cer¬ 
tain degree, slip through the horizon because the lat¬ 
ter is effectively resistive in the language of the Mem¬ 
brane Paradigm |2fi|| . in some situations the horizon 
is just not resistive enough! Indeed, the field lines are 
’’dragged” by the rotating black hole to such a de¬ 
gree that, in order for them to slip through the hori¬ 
zon steadily, they must have a certain, rather large 
toroidal field component. When, for fixed disk bound¬ 
ary conditions, the black-hole spin parameter a is in¬ 
creased beyond a certain limit Oniax, this toroidal field 
becomes so large that the poloidal field tension is no 
longer able to contain its pressure at large distances. 


IV. NUMERICAL RESULTS 

In order to verify the proposition put forward in 
the preceding section, I have performed a series of 
numerical calculations and obtained solutions of the 
force-free Grad-Shafranov equation for various values 
of two parameters: the black-hole spin a and the ra¬ 
dial extent Rg of the magnetic link. In this section I 
describe the computational set-up of the problem and 
present the main numerical results. 

The simplest axisymmetric closed-field configura¬ 
tion one could consider is that shown in Figure ^ 
where all magnetic field lines connect the disk and the 
hole. Furthermore, the entire event horizon and the 
entire disk surface participate in this magnetic link¬ 
age; in particular, the field lines threading the horizon 
very close to the axis 0 = 0 are anchored at some very 



FIG. 2: Schematic drawing of a black hole - disk magneto¬ 
sphere with a radially-limited magnetic connection. Here, 
only the inner part of the disk is coupled magnetically to 
the hole (closed field region), whereas the field lines at¬ 
tached to the outer part of the disk are open and extend 
to infinity. 


large radial distances in the disk. However, as follows 
from the arguments presented in § Hill a steady-state 
force-free configuration of this type can only exist in 
the case of a Schwarzschild black hole. This is in ac¬ 
cord with my simulations, as I was not able to obtain 
a convergent solution even for a Kerr black hole with 
the spin parameter as small as a = 0.05. 

From the discussion in §11111 we expect that, for a 
given value of a, the magnetic link between the polar 
region of the black hole and the disk cannot extend to 
distances on the disk larger than a certain i?max(a). 
Although the actual dependence i?max(a) should de¬ 
pend on the details of the problem, there are some 
general properties that it should have; in particular, 
one expects i?max to be a monotonically decreasing 
function of a and to be inversely proportional to a in 
the limit a —> 0. 

In order to test this picture and to determine the 
actual functional shape of i?max(Q), I modified the ba¬ 
sic geometry of the configuration by allowing for two 
distinct regions: the region of closed field lines con¬ 
necting the black hole to the inner part of the disk 
r < Rg, and the region of open field lines extending 
from the outer part of the disk all the way to infin¬ 
ity (see Fig. Ell. Thus, there is a critical field line, 
= '^d{Rg) < dltot, that separates open field lines 
(4> < 'I's) from the closed field lines (418 < < 41tot) 

connecting to the black hole (I count the poloidal flux 
on the disk surface from the radial infinity inward). 
Gorrespondingly, the poloidal flux on the black hole 
surface varies from 4' = It's at the pole 0 = 0 to 

= 4'tot at the equator 0 = tt/2. 

Let me now describe the computational domain and 
the boundary conditions. First, because of the as¬ 
sumed axial symmetry and the symmetry with respect 
to the equatorial plane, I performed the computations 
only in one quadrant: 0 G [0,7r/2] and r G [r//,oo]. 
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Thus, there are four boundaries: the axis 0 = 0, the 
infinity r = oo, the equator 9 = 7r/2, and the horizon 
r = rn- Of these, the axis and the equator require 
boundary conditions for 4^, whereas the horizon and 
the infinity are actually regular singular surfaces and 
so one can only impose regularity conditions on them. 

The boundary condition on the axis is simply 

4'(r, 0 = 0) = 4>s = const. (5) 

The equatorial boundary, 9 = 7r/2, actually con¬ 
sists of two parts: the infinitesimally thin disk and 
the plunging region between the inner edge of disk, 
r-m = nsco(a)j and the black hole. On the disk sur¬ 
face, r > Tin, I impose a Dirichlet-type boundary con¬ 
dition: 

'l’{r > ri^, 9 = tt/ 2) = ^dir) = ^tot —■ (6) 

r 

In the plunging region, r// < r < I set 

<r < 9 = tt/2) = 4'tot = ^'d(nn) = const. 

(7) 

This choice appears to be physically appropriate for an 
accreting disk, because the matter in this region falls 
rapidly onto the black hole and stretches the magnetic 
loops in the radial direction, greatly diluting the ver¬ 
tical field component. The horizontal magnetic field 
reverses across the plunging region, which is therefore 
represented by an infinitesimally thin non-force-free 
current sheet. 

Next, as I have discussed in § im the event horizon 
is a regular singular surface of the Grad-Shafranov 
equation and so I impose the regularity condition 0 
there. This condition is an ordinary differential equa¬ 
tion that determines the function 'I'o(0) provided that 
both and /(4>) are given. Thus, from the pro¬ 

cedural perspective, it can be viewed as a Dirichlet 
boundary condition on the horizon. 

Similarly, the spatial infinity r = oo is also a regular 
singular surface of the Grad-Shafranov equation and 
thus can also be described by a regularity condition. 
In this sense, the horizon and the infinity are equiva¬ 
lent 1^. In the particular problem set-up considered 
here, I set = 0 and /(4>) = 0 on the open 

field lines going from the disk to infinity (see below). 
Then the Grad-Shafranov equation at large distances 
is greatly simplified and the infinity regularity condi¬ 
tion yields 

'k(r = oo, 0) = dts COS0 . (8) 

In addition to boundary conditions, one has to 
specify the angular velocity of the magnetic 

field lines, which in principle should be equal to the 
Keplerian angular velocity in the disk: = 

[t’o)'!')] = •\/M/(rQ^^ -I- o'/M). However, for sim¬ 
plicity I set = 0 for the open field lines ex¬ 

tending to infinity, T < Tg. This way I don’t 



FIG. 3: The 2-D parameter space (a, ^'s). The shaded 
band running diagonally across the plot represents the 
function a^ax[Rs), the maximal value of a for which a 
force-free magnetic link can extend up to a given radial 
distance Ra on the disk. 


have to deal with the outer light cylinder that these 
lines would have to cross; correspondingly, I also set 
/('!'< dfs) = 0. In other words, I assume that the 
outer, open-field portion of the disk is not rotating 
and that the open-field part of the magnetosphere at¬ 
tached to it is potential. Next, in order to avoid the 
discontinuities in Hf('I') and /(df) at fk = ft's, I mod¬ 
ify the disk rotation law slightly just inside of Rg by 
taking Cl f smoothly to zero over a small vicinity of the 
separatrix. These modifications enabled me to focus 
on examining the effect of the black hole rotation on 
the maximal allowable radial extent Rs of the force- 
free magnetic coupling, while avoiding the numerical 
difficulties resulting from the discontinuous behavior 
of /(df), etc. I believe that these modifications do 
not lead to any qualitative change in the conclusions, 
especially in the case of small a. 

The single most important result of my calculations 
is presented in Figure El This figure shows where 
in the two-dimensional (o, ft's) parameter space force- 
free solutions exist. Filled circles on this plot repre¬ 
sent the runs in which convergence was achieved (al¬ 
lowed region), whereas open circles correspond to the 
runs that failed to converge to a suitable solution (for¬ 
bidden region). The boundary aniax(^s) between the 
allowed and forbidden regions is located somewhere 
inside the narrow hatched band that runs from the 
lower left to the upper right of the Figure. As one 
can see, aniax(^s) is a monotonically increasing func¬ 
tion and it indeed scales linearly with Ts and hence 
inversely with Rg ^ in the limit It's ^ 0, in full 

agreement with the arguments of § nni 

In order to study the effect that the black hole spin 
has on the solutions, I concentrate on several values 
of a for a fixed value of Tg. In particular, I choose 
Iks = 0.5, which corresponds to Rs = 2rin(a), and 
consider four values of a: a = 0, 0.25, 0.5, and 0.7. 


2109 



22nd Texas Symposium on Relativistic Astrophysics, Stanford, Palo Alto, December 13-17, 2004 


5 





open 


total 


FIG. 4: Contour plots of the magnetic flux function 4^(r, 9) 
for a = 0.0, 0.25, 0.5, and 0.7; the total open flux is is = 
0.54'tot in all cases, corresponding to Rb = 2rin. 


Figure 0 shows the contour plots of the poloidal mag¬ 
netic flux for these four cases. One can see that al¬ 
though the flux surfaces inflate somewhat with in¬ 
creased a, this expansion is not very dramatic, even 
in the case a = 0.7, which is very close to the crit¬ 
ical value aniax('I's = 0.5) at which a sudden loss of 
equilibrium occurs. 

After the solutions have been obtained, one can 
compute the rate of exchange of energy and angular 
momentum between the black hole and the disk. Ac¬ 
cording to M , the amount of angular momentum AL 
transported out in a unit of global time t through a 
flux tube between 'F and IF -f A'F is —/('F)A'F/2 and 
the red-shifted power is —/(iF)nir(lF)A'F/2. Then, 
since in our problem we have an explicit mapping m 
between 'F and r, we can obtain the radial distri¬ 
butions of the deposition of angular momentum and 
red-shifted energy on the disk. Figures 0 and 0 show 
these distributions for our selected cases a = 0.25, 0.5, 
and 0.7 for fixed 'F;, = 0.5. We see that in the case 
a = 0.25 there is a corotation point Tco on the disk 
such that fldisk > inside Vco and fldisk < out¬ 
side Tco. Correspondingly, both angular momentum 
and red-shifted energy flow from the inner (r < Vco) 
part of the disk to the black hole and from the hole 
to the outer (r > Tco) part of the disk. At larger 
values of a, however, the Keplerian angular velocity 
at r = Tin is smaller than the black hole’s rotation 
rate and there is no corotation point; correspondingly, 
both angular momentum and red-shifted energy flow 
from the hole to the disk. Also, as can be seen in 
Figures |S1 and ini the deposition of these quantities 



FIG. 5: Radial distribution of the magnetic torque per 
unit radius r on the disk surface for a = 0.25, 0.5, and 0.7. 



10 „ 12 


FIG. 6: Radial distribution dP/dr of the red-shifted power 
per unit radius r on the disk surface for a = 0.25, 0.5, 
and 0.7. 


becomes strongly concentrated near the disk’s edge, 
especially at higher values of a. 

A more detailed description of the numerical model 
and of the results is given in p^. 


V. CONCLUSIONS 

I investigated the structure of a force-free magneto¬ 
sphere linking a rotating Kerr black hole to its accre¬ 
tion disk. I assumed that the magnetosphere is sta¬ 
tionary, axisymmetric, and infinitely-conducting and 
that the disk is thin, ideally conducting, and Kep¬ 
lerian. My main goal was to determine under what 
conditions a force-free magnetic field can connect the 
hole directly to the disk and how the black hole rota¬ 
tion limits the radial extent of such a link. 

I first introduced a simple physical argument 
demonstrating that magnetic field lines connecting the 
polar region of a spinning black hole to arbitrarily re¬ 
mote regions of the disk cannot be in a force-free equi¬ 
librium. The basic reason for this is that these field 
lines would have to have a substantial toroidal mag¬ 
netic field component, which, in the language of the 
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Membrane Paradigm [^, is needed for the field lines 
to slip resistively across the stretched event horizon. 
In a force-free magnetosphere, toroidal flux spreads 
along the lines to keep the poloidal current constant 
along the field. Then one finds that the outward pres¬ 
sure of this toroidal field is too large to be confined by 
the poloidal field tension at large distances from the 
hole. These means that these field lines cannot be in a 
force-free equilibrium. Furthermore, one can general¬ 
ize this argument to the case of closed magnetospheres 
of finite size and derive a conjecture that the max¬ 
imal radial extent i?max of the magnetically-coupled 
region on the disk surface should scale inversely with 
the black hole spin parameter a in the limit a —> 0. 

I have then obtained numerical solutions of the 
general-relativistic force-free Grad-Shafranov equa¬ 
tion that governs the system’s behavior. To do this, 
I had to devise an iterative procedure that used 
the inner light cylinder regularity condition to deter¬ 
mine the poloidal current /(4') and the shape of the 
light cylinder simultaneously with calculating the so¬ 
lution ^'(r, 0). I performed a series of computations 
corresponding to various values of two parameters: 
the black hole spin a and the magnetic link’s radial 
extent Rg on the disk surface. I was able to chart out 
the boundary ainax('I's) between the allowed and the 
forbidden domains in the two-parameter space (a, 'ks). 
I found that this is a monotonically rising curve with 
the asymptotic behavior Umax oc d';, as 4's ^ 0, in 
line with the predictions of § nni I also computed the 
angular momentum and red-shifted energy exchanged 
magnetically between the hole and the disk. I found 


that both of these quantities grow rapidly and that 
their deposition becomes highly concentrated near the 
inner edge of the disk as the black hole spin is in¬ 
creased. 

Finally, note that, in the case of an open or 
partially-open field configuration responsible for the 
Blandford-Znajek process, one has to consider mag¬ 
netic field lines that extend from the event horizon out 
to infinity. Since these field lines are not attached to a 
heavy, infinitely-conducting disk, their angular veloc¬ 
ity cannot be explicitly prescribed; it becomes 

just as undetermined as the poloidal current /(dt) they 
carry. Fortunately, however, these field lines now have 
to cross two light cylinders, each of which being a sin¬ 
gular surface of the Grad-Shafranov equation. There¬ 
fore, one can impose corresponding regularity con¬ 
ditions on both of these two surfaces and use them 
in some coordinated, self-consistent manner to deter¬ 
mine the two free functions and simulta¬ 

neously as a part of the overall iterative scheme. 

As a final disclaimer, just in case one wonders, the 
subject of the present paper has nothing to do with 
the role of thermal conduction in galaxy clusters 
or with the excitation and propagation of eccentricity 
disturbances in planetary systems [^. 
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